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Abstract 

In order to convert the well-known Fe-Co-B alloy from a soft to a hard magnet, we propose 
tetragonal strain by interstitial boron. Density functional theory reveals that when B atoms occupy 
octahedral interstitial sites, the bcc Fe-Co lattice is strained spontaneously. Such highly distorted 
Fe-Co is predicted to reach a strong magnetocrystalline anisotropy which may compete with shape 
anisotropy. Probing this theoretical suggestion experimentally, epitaxial films are examined. A 
spontaneous strain up to 5% lattice distortion is obtained for B contents up to 4at%, which leads 
to uniaxial anisotropy constants exceeding 0.5MJ/m^. However, a further addition of B results in 
a partial amorphization, which degrades both anisotropy and magnetization. 
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I. INTRODUCTION 


Fe-Co exhibits one of the highest magnetizations among all magnetic materialsi*^ and 
is thus very attractive for applications^. Fe-Co-B alloys are well known as soft magnetic 
materials due to the glass forming ability of boron^^-. In such amorphous materials the 
magnetocrystalline anisotropy is suppressed. In this study, we demonstrate that Fe-Co-B 
is also a promising hard magnetic material, when the amorphization is avoided and the B 
atoms induce a spontaneous strain in the crystalline lattice. 

Within this introduction, we will hrst summarize criteria for soft and hard magnetic 
properties. We then use Density Functional Theory (DFT) to show that introducing B on 
interstitials sites in the Fe-Co lattice results in spontaneously strained phases with substan¬ 
tial magnetocrystalline anisotropy (MCA). In order to probe the limits of this approach 
experimentally, we study a series of epitaxial thin films. Results are discussed with respect 
to the concepts of epitaxial and spontaneous strain in Fe-Co for obtaining rare-earth free 
permanent magnets. In particular, we relate Fe-Co-B presented here to recently published 
Fe-Co-C films^. 

Since both, soft and hard magnets require a high magnetization, Fe-Co^*^ is in focus of 
research. For soft magnetic materials, the magnetic anisotropy should be as low as possi¬ 
ble, which is commonly achieved by reducing the grain size towards nanocrystalline or even 
amorphous materials^‘^‘^. Aiming at hard magnetic properties, a crystalline structure is es¬ 
sential as a magnetocrystalline anisotropy can favor particular directions of magnetization. 
The bcc crystal structure of binary Fe-Co, however, exhibits only a low cubic anisotropyl^. 
Reducing the crystal symmetry by uniaxial strain is a route to achieve high magnetocrys- 
talline anisotropy in Fe-Co as proposed in different theoretical studies^i^— . According to 
these studies, tetragonally strained Fe-Co is considered a possible rare-earth free permanent 
magnet material. Due to its high magnetostriction^'^, Fe-Co indeed appears to be more 
susceptible to strain compared to e. g. Fe-Ni, which makes Fe-Co less favorable regarding soft 
magnet applications. Though inverse magnetostriction only describes the influence of low 
strains, it allows a speculation that large strains are benehcial for a high magnetocrystalline 
anisotropy. 

The tetragonal strain is commonly described by the lattice parameters of the strained axis 
c and the compressed axes a, perpendicular to c. There are two main routes to strain the Fe- 
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Co lattice experimentally. First, in coherent epitaxial growth, Fe-Co thin hlms adapt the in¬ 
plane lattice parameter a from the substrate, which results in an induced strain (c/a > 1) due 
to compressive stress within the hlm’s plane. However, the driving forces for strain relaxation 
are too high to maintain tetragonal distortion in hlms thicker than 15 monolayersl^^— , which 
is equivalent to 2nm. In situ studies^ revealed that binary Fe-Co hlms of 4nm thickness 
are again cubic with c/a = 1. In order to stabilize the strain, calculations based on Density 
Functional Theory (DFT)^ motivated an introduction of C atoms on interstitial sites in 
Fe-Co. Delczeg et ah— proposed that about 6 at% of these small atoms establish a c/a ratio 
of 1.12. In experiments, however, the solubility on interstitial sites is much more limited. 
Using non-equilibrium preparation methods like Pulsed Laser Deposition (PLD), the limit 
of solubility, which is only about 0.1 at% for C or B in bcc Fe-Co in thermal equilibriunt^i, 
can be shifted to signihcantly higher values. PLD prepared Fe-Co hlms containing 2at% 
C exhibit such spontaneous strain with c/a = 1.03 up to at least 100 nm thickness as 
shown in our previous study^. In contrast to the induced tetragonal strain at a coherent 
epitaxial interface, interstitials can result in a minimum of total energy at a certain lattice 
distortion. This second approach is thus not limited to ultrathin hlms and accordingly, 
the strain is also present in thicker hlms. In analogy to martensitic transformations, we 
assign this distortion as spontaneous strain, though a direct examination of an associated 
transformation is difficult, as common DFT only gives the ground state and in thin hlms, a 
transformation is often hindered by the rigid substrate^. Nevertheless, the energy diherences 
calculated between cubic and tetragonal state are within a factor of two comparable to 
thermal energies at room temperature and thus transformations are expected to occur in 
bulk materials. 

Computational modeling can act as a powerful guide to experimental studies. Hence, 
we begin our work by utilizing a combination of the most relevant and highly accurate, 
spin-polarized DFT methods. In the hrst step, structural properties are modeled with low 
dopant concentrations being reproduced by using large supercells and alloying treated in 
the coherent potential approximation (CPA). The result of these structural studies are then 
used as input for the next step, where magnetic properties, including magnetic moments 
of both Fe and Co atoms on each inequivalent site, as well as MAE and orbital moments, 
are evaluated in fully relativistic calculations including the ehect of the crucial spin-orbit 
coupling. 
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In PLD, ions with high kinetic energy around 100 eV are deposited^^. Thus, besides 
regular deposition, an implantation of the material within the topmost him monolayer takes 
place^^ and allows for a supersaturation of Fe-Co hlms with C or B. However, this high 
energy impact also supports the formation of lattice defects. Supersaturated PLD prepared 
Fe-Co hlms with 2at% C retain c/a ratios of 1.03 up to high him thicknesses^. In such 
spontaneously strained Fe-Co-C, magnetocrystalline anisotropy energies (MAE) of about 
0.4MJ/m^ were observed independently from him thickness. Fe-Co with spontaneous strain 
should not have limitations in him thickness and is thus considered more promising than 
Fe-Co hlms with induced strain. 

An open question, though, is how the spontaneous strain and the MCA may be further 
increased. Since the solubility of C in Fe-Co is limited also in PLD prepared hlms^, we 
focus on boron as atom of comparable size. The location of B in bcc Fe (and similar in bcc 
Fe-Co) is still under debate. Boron can substitute Fe in the ideal bcc positions^l^iSS^ can 
occupy interstitial positions^^— , and B solution in Fe can be stabilized by Ni, N and C^i^ 
both, interstitially and substitutionally^. In a more recent theoretical investigation^, it was 
shown that B prefers the substitutional chemical disorder in bulk Fe-B. However, this work 
did not take a local relaxation of the strain around interstitial B atoms into account. The 
authors further found that surfaces stabilize the occupation of octahedral interstitials. Near 
free surfaces, the geometrical pressure on B interstitials relaxes. Occupation of octahedral 
interstitials is followed by an increase of the B-Fe distances and thus a distortion of the 
bcc lattice. Baik et ah— thus give the motivation to further study the straining effect of B 
interstitials in Fe(-Co) thin hlms which are also dominated by free surfaces. 

Considering the different interstitial sites in a bcc lattice, a tetragonal strain is exclusively 
expected, when the B atoms occupy octahedral interstitials. Substitutional arrangements 
or occupation of tetrahedral interstitials would not change c/a, but only affect the unit cell 
volume. Octahedral interstitials as dehned in Ref.— exist along all three spatial directions 
(see supplementary^^). There are six octahedral interstitial sites in a bcc lattice, indicating 
a high theoretical solubility, when the atomic radius of the interstitial atom is sufficiently 
low and chemical conditions are benehcial. If such an interstitial site is occupied by a small 
atom like B, the atoms of the apexes of the octahedron displace along its axis, i. e. uniaxially. 
A lattice strain by means of c/a > 1 is thus only possible, when the octahedral interstitials 
along the c axis as e.g. (0;0;l/2) within the bcc unit cell are preferentially occupied. The 
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question, how an additional occupation of the other octahedral institials affects the structural 
properties and the total energy will be discussed based on DFT calculations. 

Besides the discussed lattice preferences of boron, amorphization of the Fe-Co lattice is 
expected to begin at a certain B content. Depending on the film preparation conditions, 
amorphous Fe-Co-B phases were reported at B concentrations of e.g. 7.5^, 15^ or 22at%^. 
In order to establish a high magnetocrystalline anisotropy in Fe-Co-B, amorphization has to 
be avoided. Our study thus presents theoretical calculations of the preference of B atoms to 
occupy the same type of octahedral interstitials due to lower energy of such conhguration. 
Then it focuses on the MCA evaluation of ideal Fe-Co-B crystals assuming all B impurities 
occupying the octahedral interstitials perpendicular to the substrate. We then introduce the 
properties of PLD prepared epitaxial Fe-Co-B films, taking the different possibilities how B 
may affect the Fe-Co lattice into consideration. 


II. METHODS 

A. Computational methods 

First principles electronic structure calculations within the generalized gradient approxi¬ 
mation (GGA)^ were used to identify stable or metastable body centered tetragonal (bet) 
structures for B doped Fe-Co alloys. 

In the first step, different B concentrations were modeled by three Fe^B supercells, FegB, 
FeieB and Fe 24 B with 1x2x2,2x2x2 and 2x2x3 supercells, respectively. Fe atoms 
occupy the ideal bcc positions, while B atoms are placed in the octahedral positions, which 
are the most stable interstitial positions^. B dopants in the octahedral positions reproduce 
the upper limit for tetragonal distortion. Substitutional B dopants were not considered as 
they do not contribute to the formation of distortion. The conjugate-gradient algorithm 
as implemented in the Vienna Ab Initio Simulation Package (VASP)^^— was used to fully 
relax these structures. The k-point mesh was set to 16 x 16 x 8, 8 x 8 x 8, 8 x 8 x 6 for 
FegB, FeigB and Fe 24 B, respectively, within the Monkhorst-Pack scheme^. The plane-wave 
cut off and the energy convergence criterion of the scalar relativistic calculations were set 
to 500 eV and 10“^ eV, respectively. 

In the second step, the alloying effect on the equilibrium parameters (c/a)eq-ratio and 
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volume Veq of these relaxed structures was studied using the Exact Muffin-Tin Orbitals - 
Full Charge Density (EMTO-FCD) method^^— . The accuracy of the EMTO method is 
controlled by the optimized overlapping muffin-tin potential^‘^>^. The muffin-tin potential 
optimization procedure is described in details in Ref.—. The potential parameter rj was cho¬ 
sen to match the equilibrium parameters obtained by EMTO to those obtained by the VAST 
calculations. The potential parameter rj is 0.75, 0.85 and 0.9 for FegB, FeieB and Fe 24 B, 
respectively. The chemical disorder for (Fei_ 3 ;Coj;)yB alloys was treated via the Coherent 
Potential Approximation (CPA)^!^. The one-electron equations were solved within the 
scalar-relativistic and soft-core approximations. The Green’s functions were calculated for 
18 complex energy points distributed exponentially on a semi-circular contour with radius 
of 1.2 Ry. The 3d and 4s states of Fe, Co, and the 2s and 2p states of B were treated as 
valence electrons, s, p, d, f orbitals were included in the muffin-tin basis set (/max = 3). 
The one center expansion of the full-charge density was truncated at = 8. Between 40 
and 500 uniformly distributed k-points were used in the irreducible wedge of the Brillouin 
zone. The electrostatic correction to the single-site CPA was described using the screened 
impurity model^*^ with a screening parameter of 0.6. Other screening parameters^ have 
been tested for the smallest system and it turned out that their effect on the equilibrium 
parameters is less than 1 %. 

The magnetic properties, including magnetic moments and the magnetocrystalline 
anisotropy energy (MAE), of the various (Fei_ 2 ;Co 3 ;)j^B systems were evaluated using the 
spin-polarized relativistic KKR (SPR-KKR)^!^ method in the atomic spheres approxi¬ 
mation (ASA), in a similar manner as was done in Ref.—. Calculations were performed 
using GCA— for the exchange-correlation potential and at least several thousand k-points 
(depending on system size) were used for numerical integration over the Brillouin zone in 
order to obtain well converged values of MAE. Alloying was treated with the GPA— and 
the MAE was obtained by total energy difference for two different magnetization directions, 
i. e. MAE = Etotai(ffi || 100) - Etotai(ffi || 001). 

All structures from subsection IIIIBI were fully relaxed, c/a ratios and volumes were 
optimized together with atomic relaxation for every c/a ratio and volume. Those calculations 
were performed with the full potential linearized augmented plane wave method (FP-LAPW) 
implemented in the WIEN2k code^. For the exchange-correlation potential the GGA form^ 
was used. Galculations were performed with a plane wave cut-off parameter RiCmax = 6.5, 
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total energy convergence criterion 10“® Ry, with 4x4x4 k-points and with radii of the 
atomic spheres 1.96 oq for Fe/FeCo and 1.45 oq for B, where oq is the Bohr radius. The 
Virtual Crystal Approximation (VGA) was used to study different Feo. 4 Coo. 6 -B compositions, 
following the same procedure as described in Ref.—. All parameters were carefully tested to 
provide well converged values. 


B. Experimental 

The Fe-Co-B samples were prepared as thin dims performing Pulsed Laser Deposition 
(PLD) in ultra-high vacuum (5 x 10“®mbar) at room temperature. The used KrF excimer 
laser (Coherent LPXpro 305) operates with a wave length of 248 nm and a pulse length of 
25 ns. On MgO(lOO) single crystal substrates, 3nm Cr seed layers and 30 nm Auo.55Cuo.45 
buffer layers were deposited prior to the deposition of the 20 nm thick Fe-Co-B dims. The 
Au-Cu and the Fe-Co-B layers were prepared in pseudo-co-deposition, i. e. by repetitive 
changing of the targets during PLD to achieve the aimed compositions. Elemental (An, 
Cu, Fe and Co) and a FeB composite target were used therefore. The deposition rates were 
measured prior the preparation with a quartz crystal rate monitor. 

Energy dispersive x-ray spectroscopy (EDX) measurements on a Bruker EDX in a JEOL 
JSM6510-NX electron microscope condrmed the intended dim compositions. The boron 
content in the Fe-Co-B dims was measured exploiting Auger Electron Spectroscopy (AES) 
on a JEOL JAMP-9500F Field Emission Auger Microprobe device. The atomic concentra¬ 
tions^ were calculated using standard single element sensitivity factors from mean values 
of sputter depth prodles measurements carried out with 1 keV Ar"*" ions. Atomic force mi¬ 
croscopy (AFM) measurements on an Asylum Research Cypher AFM were performed to 
study the surface morphologies and to condrm a high datness of the dim surfaces. A Bruker 
D 8 Advance didractometer operating with CoKa radiation was used for X-ray didraction 
(XRD) in Bragg-Brentano geometry. The Fe-Co-B lattice strains were derived from pole dg- 
ure measurements, which had been carried out on an X’pert four circle goniometer (CuKa 
radiation). Transmission electron microscopy (TEM) studies were performed on a Titan^ 
80-300 microscope, which was equipped with a Cs corrector and a Schottky held emission 
electron source. The lamella was cut using focused ion beam on a FEI microscope (Helios 
NanoLab 600i). For magnetic hysteresis measurements at 300 K, a vibrational sample mag- 
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netometer (VSM) mounted in a Quantum Design Physical Property Measurement System 
was used. 

III. RESULTS AND DISCUSSION OF DFT STUDIES 

A. Structure optimization 

In order to prepare crystal structures of chemically disordered (Fea:Coi_ 2 ;)yB, the simpler 
Fe^B systems have to be considered as an initial step. Fe^^B structures are modeled with 
B atoms placed into octahedral interstitial positions along the c axis. This configuration 
determines the upper limit of the tetragonal distortion. Internal relaxation of Fe atoms 
around the B interstitials leads to the elongation of c and shrinking of a lattice parameter. 
In other words, the c/a ratio is now higher than 1, where 1 corresponds to a cubic lattice. 
The total energies E as a function of c/a ratio and volume V for three Fej^B systems are 
computed by VASP and their plots are presented in Figured] 


FIG. 1. (Color online) Contour plots of total energy (meV/Fe atom) as a function of c/a-ratio and 
unit volume V, for FegB, FeigB and Fe 24 B, respectively. The zero reference values for the energies 
are E{c/a=l, The zero energy isolines are denoted by dashed lines. Total energy and 

volume are normalized per Fe atom. 


At each c/a ratio, the equilibrium volume (I4q) is obtained by using a Morse-type equation 
of state as minimum of E{V) dependency. Next, the equilibrium ratio (c/a)eq is attained 
by htting E{Veq) with a second order polynomial. Results of (c/a)eq and Ueq are listed 
in Table E It is observed that both the c/a ratio and volume increase with B content, as 
expected. 

Three types of (Fea;Coi_a;)yB supercells corresponding to different B concentrations have 
been considered, i. e. (Fe 3 ;Coi_ 2 ;) 8 B, (Fea;Coi_^)i 6 B and (Fe^Coi_ 3 ;) 24 B, see Figured] Ini¬ 
tially, (Fe 3 ;Coi_a;) 8 B and (Fea;Coi_ 3 ;)i 6 B were taken into account as reasonable crossing point 
between computational effort and expected experimental B solubility. Preliminary results 
suggested that higher B concentration gives higher c/a ratio and MAE. Unfortunately, nei¬ 
ther 11.1 at% B ((Fe 3 ,Coi_ 3 ;) 8 B) nor 5.9 at% B ((Fe 3 ;Coi_ 3 ,)i 6 B) could be obtained experi- 


TABLE I. Equilibrium (c/o)eq ratios and volumes (14q) of Fe^^B systems obtained by EMTO and 
VASP. 


Alloy 

at% B 

(c/a)eq 

VASP EMTO A 

I 4 q(A^/Fe atom) 

VASP EMTO 

FegB 

11.1 

1.126 1.132 -0.006 

12.73 12.88 

FeieB 

5.9 

1.089 1.084 0.005 

12.21 12.25 

Fe24B 

4.0 

1.050 1.044 0.006 

11.93 11.95 



FIG. 2. (Color online) Crystallographic models of (Fe 3 ;Coi_ 3 ;) 8 B, (Fe 2 ;Coi_ 3 ;)i 6 B and 
(Fe 3 ;Coi_a;) 24 B, respectively. The bigger red spheres represent Fe/Co atoms while smaller black 
spheres denote B atoms. 


mentally, and what was shown instead, maximum B solubility producing maximum tetrag¬ 
onal strain is about 4at%. At that point of this study the (Fea;Coi_ 3 ;) 24 B structures (4at% 
B) were introduced and calculated. 

The effect of substituting Fe by Co on the structural parameters of Fe^B has been taken 
into account by using the CPA approximation implemented in the EMTO method. The 
equilibrium ratios (c/a)eq and the lattice parameter Oeq have been evaluated for the rigid 
supercells taken from the preceding VASP calculations. Thus, it is important to assess the 
accuracy of the EMTO method compared to VASP for the pure Fe^B. For each c/a ratio, 
a fully relaxed structure is generated by VASP at lAq and used as input for the EMTO. 
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Total energies are then calculated for each c/ a ratio and five different V's while keeping the 
internal parameters fixed. Results obtained by EMTO, using the above mentioned fitting 
procedure for VAST, are listed in Table HI The error bar of the (c/a)eq ratio is defined by A 
= (c/a)y^®^ - (c/a)™'^°, with c/a ratios obtained by VAST and EMTO. Good agreement 
between these two theoretical methods, regarding both (c/a)eq and 14 q, is observed. 

After validation of the EMTO method, the desired CPA structures of (Fea;Coi_ 3 ;)j,B have 
been evaluated. Results are presented in Table HIl R is observed that an increase of Co 
concentration in (Fea;Coi_ 3 ;)j,B enhances the tetragonal distortion (c/a)eq. At the same time 
the Oeq and Kq decrease. Decrease in Veq is ascribed to the lower atomic volume of Co 
compared to Fe. 

Note, that the tetragonal distortion predicted for B-doped Fe-Co alloys is larger than for 
C-doped alloys^^. This can be understood by the fact that the B atomic radius is bigger 
than the C atomic radius. By adding 4at% of B a distortion around 1.05-1.07 should be 
achieved, while for C-doping it is merely about 1.03-1.04. This comparison suggests that 
B-doping could be a better way to achieve tetragonally distorted Fe-Co alloys. 

Note further, that with B-doping, the tetragonal distortion can be obtained for a broader 
range of Co concentrations, compared to C-doped counterparts. This suggests a possibility 
to prepare tetragonal Fe-Co alloys with lower Co contents. In order to find out whether this 
brings any advantages, we discuss their magnetic characteristics in the following section. 


B. Preferential orientation of octahedral interstitials 

For a permanent magnet it is decisive that the easy axis of all unit cells is aligned along 
one particular direction, in which the magnet will be magnetized. The different orientations 
of octahedral interstitials^^, which could be occupied by the B atoms, thus have to be consid¬ 
ered. Consequences of a random distribution have to be discussed as well. The occupation 
of tetrahedral sites, which just result in a volume change, are not considered in the following 
since they are energetically unfavorable. Baik et ah— computed that the energy difference 
between octahedral and tetrahedral sites amounts to 0.70 eV for a cubic Fe-B supercell, with¬ 
out taking into account an optimization of lattice parameters. The corresponding energy 
difference calculated by us is 0.77 eV/B atom, based on bcc Fe 2x2x2 supercells (FeieB) with 
geometry optimization. This means, that B atoms favor the octahedral sites even stronger 
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TABLE II. Summary of calculated structural (EMTO) and magnetic (SPR-KKR) characteristics 
of (Fe 3 ;Coi a;)jyB systems. 


Composition 

(c/a)eq 

Cteq 

(A) 

Feq 

( i 

ms 

( ) 

mi 

i i 

Ho Ms 

(T) 

MAE 

/ fj.eV \ 

MAE 

tMJt 

'' Fe atom 

atom ^ 

'' atom ^ 

^ atom ^ 

(Feo.5Coo.5)8B 

1.247 

2.716 

12.49 

1.71 

0.066 

1.87 

118 

1.51 

(Feo.5oCoo,5o)i6B 

1.086 

2.809 

12.04 

1.97 

0.068 

2.09 

46 

0.62 

(Feo.45Coo,55)l6B 

1.091 

2.802 

12.00 

1.91 

0.068 

2.04 

51 

0.69 

(Feo.4oCoo,6o)i6B 

1.103 

2.788 

11.95 

1.86 

0.069 

2.00 

52 

0.69 

(Feo.35Coo.65)l6B 

1.116 

2.774 

11.91 

1.82 

0.069 

1.96 

43 

0.58 

(Feo.esCoo, 35)246 

1.049 

2.835 

11.95 

2.16 

0.064 

2.26 

23 

0.31 

(Feo.eoCoo. 40)246 

1.051 

2.832 

11.94 

2.12 

0.066 

2.22 

30 

0.40 

( 600 . 55000 . 45)246 

1.053 

2.827 

11.90 

2.08 

0.067 

2.19 

30 

0.40 

(Feo.5oCoo.5o)246 

1.056 

2.822 

11.87 

2.04 

0.068 

2.15 

28 

0.38 

( 600 . 45000 . 55)246 

1.059 

2.816 

11.82 

1.99 

0.069 

2.12 

27 

0.37 

(60o.4oCoo.6o)246 

1.063 

2.810 

11.79 

1.95 

0.070 

2.08 

29 

0.40 

( 600 . 35000 . 05)246 

1.068 

2.803 

11.76 

1.91 

0.070 

2.04 

35 

0.48 


than Baik et al. showed. A similar set of calculations for the alloy (Feo. 4 Coo. 6 )i 6 B, treated 
within the VGA, gives the energy difference 0.39 eV/B atom, which also suggests to consider 
only the octahedral interstitials in the subsequent models. 

As our calculations reveal a tetragonal distortion with uniaxial anisotropy along the c 
axis, independent of the B content, we can take any tetragonal building block for the follow¬ 
ing considerations, which are borrowed from martensite theory as this can be extrapolated 
towards the unit cell level^*^. We consider martensite theory as reasonable since the diffu¬ 
sion of B as interstitial^ is much faster in contrast to Fe (or Co) atoms^. In this picture, the 
orientation of the tetragonal distortion of the Fe-Co ’’martensite” changes, when B moves 
from one octahedral site to another. 

The connection of different tetragonal building blocks is possible by twin boundaries, but 
they require additional twin boundary energy. A twin boundary is connected with a 90° do¬ 
main wall in these uniaxial ferromagnets, which further requires magnetic exchange energy. 
In this simplihed picture, both, elastic and magnetic contributions make a different orienta- 
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tion of neighboring tetragonal building blocks energetically unfavorable. Twin boundaries 
must be introduced, however, to minimize elastic stress energy and magnetostatic energy, 
and in bulk samples all three orientations are thus equivalent and should occur together. 



a) Fe 54 B 2 orthogonal b) Fe 54 B 2 uniaxial 


FIG. 3. 3x3x3 bcc Fe supercells with two B atoms in octahedral interstitial positions (Fe 54 B 2 ) as 
a) orthogonal and b) uniaxial configuration of interstitial axes. 

In order to support this picture on the atomic scale, we performed DFT calculations of 
a conhguration with orthogonal B interstitials and compare it with a uniaxial orientation of 
B interstitials. For this comparison, one has to consider at least two B atoms per unit cell, 
whereby the Fe/FeCo supercell has to be enlarged in order to keep low B concentration. In 
the hrst step, two Fe based models are considered, both with two B atoms in octahedral 
interstitial positions in bcc Fe 3x3x3 supercells (Fe 54 B 2 ), see Fig. |3l In the orthogonal 
model (Fig.|3]^a)), the positive strains are generated along two orthogonal axes, which leads 
to a c/a ratio below 1. In the second model (Fig.|3]|^b)), the octahedral interstitial axes 
are oriented parallel, which forms strain exclusively along the c axis and leads to c/a ratio 
above 1. Calculations show that the uniaxial conhguration is favored with a difference 
of total energy equal to 129 meV/B atom. In order to relate the theoretical result to the 
experimental Fe-Co-B system one has to go beyond the pure Fe model and, as we did in the 
previous section, consider the inhuence of alloying with Co. For this reason, in the second 
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step, the two Fe 54 B 2 models described above have been recalculated for a selected Feo. 4 Coo .6 
concentration within the VGA. Results show once again that the uniaxial conhguration is 
more stable, with slightly lower total energy difference of 126 meV/B atom. CPA calculations 
are prohibitively computationally expensive for the EMTO method due to a large number of 
inequivalent atoms in the unit cells. Some crystallographic data of the optimized structures 
are collected in Table m An interesting result is that the supercell volumes, which for 
orthogonal and uniaxial conhgurations were optimized independently, have the same values 
up to 5 signihcant digits, which indicates the high accuracy of these results. 

TABLE III. Crystallographic data of optimized Fe 54 B 2 and (Feo. 4 Coo. 6 ) 54 B 2 structures and differ¬ 
ences of total energies between orthogonal and uniaxial configurations of B interstitials. 



orthogonal 

uniaxial 



V a c/a 

V a c/a 

AE 


[A3] [A] 

[A3] [A] 

meV 


B atom 

Fe54B2 

635.74 8.652 0.982 

635.74 8.492 1.038 

129 

(Feo,4Coo.6)54B2 

625.31 8.636 0.971 

625.31 8.454 1.035 

126 


With three B atoms per 3x3x3 bcc Fe supercell it is possible to construct a system with 
octahedral interstitials in each of the three spatial directions, resulting in a cubic structure. 
Such a Fe 54 B 3 structure with a simulated random occupation of the octahedral interstitials, 
while keeping their mutual distances as large as possible within the supercell, has been 
carefully optimized together with its uniaxial Fe 54 B 3 counterpart. The total energy differ¬ 
ence between these structures is 49meV/Batom and the preferred configuration of boron 
is again the uniaxial one. After dividing the latter value by Boltzmann constant one ob¬ 
tains 570 K/B atom, the same order of magnitude as room temperature. This implies, that 
a small fraction of B atoms will occupy positions not yielding a tetragonal strain. Never¬ 
theless, calculations of the presented ideal cases, which may not compare the energies of 
all theoretical possibilities, but compare the most extreme scenarios, suggest that an occu¬ 
pation of neighboring orthogonal octahedral sites with varying orientation is energetically 
unfavorable. 

From an experimental point of view, the preferential site occupation is tedious to be 
measured directly due to the much lower atomic number of B compared to Fe-Co. However, 
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the reduced symmetry of thin dims compared to bulk may result in a favorable orientation, 
which could be induced by coherent epitaxial growth, followed by a relaxation towards the 
spontaneously strained state as shown for Fe-Co-C films^. A preferential lattice orientation 
allows using integral methods like texture measurements to probe the present tetragonal 
distortion and magnetization measurements to determine the magnetocrystalline anisotropy. 
An agreement of these global measurements with the local DFT calculations would also 
conhrm the unique alignment of all easy axes, which is benehcial for permanent magnet 
applications. 

In other words, when comparing DFT calculations with experiments, one has to consider 
the different length scales. In DFT calculations, the preferential c axis is the result of 
the spontaneous strain at the atomic scale, while in thin film experiments, the out-of-plane 
orientation is given and may act on the whole sample. The expected preferential alignment of 
the strained c axis along the out-of-plane orientation thus has to be probed by the following 
experiments. 


C. Saturation magnetization 

For the considered (Fea;Coi_a;)j,B systems, the influence of B on magnetic moments is re¬ 
garded first. The average magnetic moments per atom, obtained by SPR-KKR calculations, 
are listed in Table HIl R is observed that the average magnetic moment decreases by increas¬ 
ing the B content. This can be simply explained by addition of a non-magnetic B component 
itself, and furthermore by the reduction of the magnetic moments on Fe/Co atoms around 
the B interstitial. In order to look closer at the latter effect an insight into atom and site 
specihc magnetic moments is necessary. Internal relaxation changes the distances between 
Fe/Co atoms, especially around B interstitials. These reconhgurations affect the exchange 
coupling between Fe/Co atoms and thus affects magnetic moments. 

In Table HVl the local magnetic moments for selected (Feo. 5 Coo. 5 )j,B systems are presented. 
The Fe/Co atoms that are the closest neighbors to the B atoms have the magnetic moments 
most rednced compared to the moments on Fe/Co atoms far from B. The latter values alter 
rarely, indicating that the Fe/Co atoms far from B have bulk-like characteristics already for 
(Feo. 5 Coo. 5 ) 8 B. This suggests that the influence of a B impurity is relatively short-ranged. 

An increase of Co content leads to a decrease of the total magnetic moment for the 
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TABLE IV. Magnetic moments of (Feo, 5 Coo. 5 )j/B systems (/XB/atom) obtained by SPR-KKR. Mo¬ 
ments on Fe and Co atoms in the nearest neighborhood of B octahedral interstitial and as far as 
possible from B. 


Alloy 


close to B 

far from B 


at% B 

m-Fe mco 

mpe m-co 

(Feo,5Coo.5)8B 

11.1 

1.74 0.98 

2.89 1.95 

(Feo.5Coo.5)i6B 

5.9 

1.52 0.82 

2.85 1.93 

(Feo.5Coo,5)24B 

4.0 

1.52 0.85 

2.76 1.91 


considered Co concentrations in the (Fea;Coi_ 2 ;) 24 B system (see Table m- This negative 
dependence is not necessarily true beyond the studied range, since for bet Fe/Co alloys a 
magnetization maximum is observed for lower Co concentrationsi>^. In our case the drop of 
magnetic moment is accompanied by an increase of (c/a)eq. 


D. Magnetocrystalline anisotropy energy 

Table UTI summarizes the MAEs of the (Fe 3 ;Coi_ 3 ;)yB systems, as obtained by SPR-KKR 
calculations. Trends are in many aspects similar as those obtained for interstitial carbon^. 
In particular, the tetragonal distortion of the crystal structure leads to signiheant values 
of the MAE under particular c/a and alloy concentrations. Here, however, the reduction 
in MAE is rather small when going from y = 16 to y = 24, most likely - as pointed out 
above - because the c/a values are larger for y = 24 with B than with C, which leads to 
larger MAEs compared to similar dopant concentrations of C. Considering the MAE as a 
function of Fe concentration x for the case oi y = 24, there is a rather flat behavior, but 
a maximum value of 0.48 MJ/m^ is obtained for x = 0.35. For y = 24 and x = 0.40 we 
obtain MAE = 0.40 MJ/m^, which can be compared to MAE = 0.19 MJ/m^ obtained with 
C-doping instead of B. This means that a rather small amount of interstitial B dopants 
appears to be enough to yield a signiheant MAE which is promising in a permanent magnet 
context. 

For systems with y = 24, we have performed test calculations (not shown) utilizing a 
simpler approach to the problem of alloying, the virtual crystal approximation (VCA), fol- 
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lowing the same procedure as described in Ref.—. Although these calculations overestimate 
the MAE by a factor between 2 and 4 compared to CPA calculations, nevertheless they 
result in a similar trend and thus provide a further support to the Ending that a substantial 
MAE can be obtained also for lower Co contents. In contrast to CPA results, VCA ones 
were obtained within a full potential method. 


IV. EXPERIMENTAL RESULTS AND DISCUSSION 

A. Structural properties of Fe-Co-B films 

Using PLD, we prepared a boron composition series and studied the structural and 
magnetic properties of Fe-Co-B films with 20 nm thickness. All films were deposited on 
Auo. 55 Cuo ,45 buffer layers, which provide a reduced in-plane lattice parameter compared to 
Fe-Co and thus induce a c axis oriented film growtb^>^. Based on the theoretical results, 
where the Fe to Co ratio was varied, a composition close to Feo. 4 Coo .6 was chosen as base for 
our experimental investigations. For the films of this (Fe 3 ;Coi_ 3 ;)-B series, EDX measure¬ 
ments gave X = 0.38(2). Corresponding XRD patterns are shown in Figure 111(a). The (002) 
reflection is the only intensity originating from the Fe-Co-B films and indicates an epitaxial 
film growth of Fe-Co-B on Au-Cu. No additional phases as e. g. borides, are detected. Pole 
figure measurements (see supplementary^^) confirm the epitaxial growth without formation 
of twins. Compared to the position of the Feo. 38 Coo .62 (002) reflection given in literature^, 
we observe a shift to lower 26 angles for all samples. This indicates already an expanded c 
axis of the Fe-Co-B lattice, which will be discussed later in detail. 

The AES depth profiles, which were carried out to determine the B content of the Fe-Co 
films, did not reveal a significant variation of the B content in dependence on the sputter 
depth. This allowed us to determine error bars of the B contents for each considered film from 
its statistical variation. With increasing B content, both a reduced scattering intensity and 
a peak broadening of the Fe-Co-B(002) reflection is detected. Applying Scherrer’s formula, 
we estimate the x-ray coherence length of the Fe-Co-B films, which can be taken as measure 
for the Fe-Co-B crystal size. The results are plotted in Figure ll](b). An increase in B content 
leads to a substantial decrease of coherence length from almost 10 nm in binary Feo. 38 Coo .62 
to about 3nm for the film with 9.6 at% B. 
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FIG. 4. (Color online) (a) XRD patterns of 20 nm Feo, 38 Coo. 62 -B films with varying B content on 
Auo. 55 Cuo .45 buffers. Relevant reflections are marked. The ( 002 ) equilibrium 29 angle of binary 
Feo. 38 Coo .62 is given as broken line. The XRD pattern of a Feo. 38 Coo .62 film without B addition is 
added for comparison. The intensity maxima near 60° and 77° are attributed to PLD droplets with 
the composition of the used targets, (b) shows the x-ray coherence length according to Scherrer’s 
formula against the AES determined B content. 

A TEM study of a Fe-Co-B film with 4.2 at% B confirms the crystallinity and continuity 
of the film (Figure IS]). The Fourier Transform (FT) of the Fe-Co-B film (inset) gives a very 
regular pattern as expected for an epitaxially grown crystal. Minor contrasts within the 
Fe-Co-B film may indicate the formation of separated grains with a grain size of approx. 
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FIG. 5. (Color online) TEM image of a 20 nm Feo. 38 Coo .62 film with 4.2 at% B deposited on 
Auo. 55 Cuo. 45 . The zone axis is [010]. The inset shows a FT image of the Fe-Co-B layer. The 
marked reflections were used for determination of the lattice parameters c and a. 

5nm. This would be in agreement with the observed x-ray coherence length (Figure |l](b)), 
but may also be linked to thickness variations of the lamella. 

The strong reduction of crystal size at higher B contents as indicated by the XRD mea¬ 
surements (Fig. 111(b)) has been confirmed by an additional TEM investigation of the film 
with the highest B content (see supplementary^^), which revealed oriented nanocrystals of 
the same diameter surrounded by nanocrystalline and amorphous phases. In literature, 
both, a formation of Fe-Co-B nanocrystals^ or amorphous Fe-Co-B^ is reported for films of 
similar B contents. Asai et al. also observed an amorphous structure already at B contents 
above 5at% in addition to crystalline Fe-Co-B in Feo. 7 Coo. 3 -B films^. Our measurements 
thus confirm their findings of a decreasing size and fraction of Fe-Co-B crystals with in¬ 
creasing B content in the films. However, all of the presented films contain a significant 
crystalline fraction which is grown epitaxially. These crystals comprise a certain amount of 
the alloyed B atoms as is clearly indicated by the change of their lattice parameter c, when 
compared to the binary Feo. 38 Coo .62 film, which is given as reference in Figure ID) a). 

In order to characterize the tetragonal strain in the epitaxially grown Fe-Co-B films, 
{011} pole figures were measured for all films. As described in Ref.-, the c/a ratios were 
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FIG. 6 . (Color online) (a) Tetragonal distortion c/a for 20 nm thick Feo. 38 Coo. 62 -B films in depen¬ 
dence of the B content (full squares). Results of Fe-Co-C films^ (open diamonds) have been added 
for comparison. Note that the film with the highest C content is of Fe-C type and did not contain 
Co. (b) Corresponding uniaxial magnetic anisotropy constants K\] of the films. The results from 
the DFT calculations Feo. 4 Coo. 6 -B (Table HI]) have been added as open stars in both graphs. Cubic 
Fe-Co [c/a = 1) with no uniaxial MCA [K\] = 0) has been added, although it was not treated in 
our calculations. Broken lines are guides to the eye. 

determined from the tilt angles of these planes. In all dims, one predominant variant with 
longer c than in-plane axes was found. No indications for variants with deviating c axis 
orientation were found. The c/a ratios are plotted in Figure | 6 ](a) for the different boron 
contents. An increase of tetragonal distortion is observed with increasing B content up to 
4.2 at%. The maximum observed c/a ratio is 1.045 and is achieved for this composition. We 
can compare this hnding to Feo. 4 Coo. 6 -C dims (open diamonds), which exhibit a maximum 
tetragonal strain of c/a = 1.03 and a lower dependence on the C content as discussed in our 
previous study^. The observation of an increased strain in the Fe-Co-B dims qualitatively 
condrms our theoretical calculations, which predicted an enhanced c/a ratio by about 0.03 in 
Fe-Co-B when compared to Fe-Co-C. The maximum observed strain with c/a of nearly 1.05 
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is also reasonably close to the predicted c/a = 1.063 for Feo. 4 Coo .6 with 4at% B fTable [IT|) . 
However, a further increase of B content is not followed by a further increased tetragonal 
strain, which was predicted by DFT. There are different possible reasons for this difference: 
(a) The limited solubility of boron in Fe-Co which is nearly zero in thermal equilibriunt^i, but 
is obviously increased substantially by PhD, an effect known for many systems^. A limited 
solubility favors nanocrystalline or amorphous B-rich phases which reduce the effective B 
content in the Fe-Co matrix and thus also decreases the tetragonal strain, (b) The possible 
substitutional alloying of B as discussed within the introduction, which would beneht a cubic 
lattice. And/or (c) the possibility that not all B atoms occupy octahedral interstitials along 
the axis perpendicular to the him surface, but also interstitials along the in-plane axes are 
occupied, as e. g. (l/2;0;0) within the bcc unit cell, which beneht a certain in-plane strain, 
i. e. a reduction of the measured c/a. 

All reasons have in common that the amount of B atoms, which contribute in a lattice 
strain along the c axis, is ehectively reduced. Our DFT calculations were based on the 
assumption that all B atoms occupy octahedral positions along the c axis. As shown in the 
calculations in section UlIBl the B atoms indeed preferentially align in octahedral interstitials 
along one particular axis, although all three crystal axes are equivalent in a Fe-Co lattice. 
The experimental results thus conhrm these calculations. However, because of the named 
reasons, not all B atoms contribute to the c axis oriented strain above a B content of 4.2 at%. 
A likely explanation, why DFT calculations are not able to conhrm the experiments exactly 
is, that they do not consider kinetic ehects. Films are prepared within a hnite time, which 
may lead to local variations of B density and thus result in an occupation of diherent lattice 
sites than expected from DFT, which describes the ground state. 

The comparison of the B and C doped Fe-Co hlms^ implies that a higher amount of B 
than C atoms can be solved in the PhD prepared Fe-Co hlms and preferentially occupies 
interstitial sites along the c axis, which favors a higher strain. An additional contribution 
may come from the size of the B atoms, which are bigger than C atoms. However, when a 
certain limit is reached, further added B (or C) atoms lead to a decrease of c/a. This limit is 
reached at around 4at% B (or C). Since the c axis length is not altered when adding more B 
or C (see Figure 0]), we conclude that most of the additional atoms occupy sites along the a 
axes or are dissolved from the Fe-Co crystals. The latter explains the decreasing crystal size 
with increasing B (or C) content as discussed with Figure |l](b) and in Ref.-, respectively. 
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For the Feo. 38 Coo .62 him with 4.2 at% B, the lattice distortion was determined from the 
FT of the TEM images (Figure [5]). Compared to the XRD result, where c/a was 1.045 
(Figure [Sl(a)), we observe a reduced strain of 1.02±0.01. Such strain reduction after TEM 
lamella preparation was already reported for Fe-Co-C Elms and was attributed to a lattice 
relaxation due to a reduced constraint of a thin TEM lamella compared to a continuous 
hlm^. 

From growth studies of Fe-Co-C Elms, we concluded that the misEt dislocations which 
form during Elm growth afiect the whole underlying Elm^. Most of the crystal volume is 
thus expected to exhibit the same strain state. Only a small fraction of the Elm close to 
the bufier interface had an increased strairn^, expressed by a c/a ratio increased by 0.02. 
This also holds for the Fe-Co-B Elms examined here. The tetragonal strain determined 
close to the Au-Cu bufier from TEM images of the 4.2 at% B containing Feo. 38 Coo .62 Elm 
is 1.05±0.04. This is about 0.03 higher than the value for the Elm’s volume as determined 
in TEM, which may be explained by the lower in-plane lattice parameter of Auo.55Cuo.45 
compared to Fe-Co^. The TEM studies thus prove that only a very small fraction of the 
Elms is inEuenced by this induced strain at the bufier-Elm interface and most of the Elm 
exhibits a spontaneous strain. 

The structural characterisation of the Elms clearly supports our picture of tetragonally 
strained Fe-Co-B lattices, where the axis perpendicular to the Elm surface is preferentially 
strained and the in-plane axes are compressed. The decisive question, why a predominant 
amount of B atoms preferentially occupies the octahedral interstitials along this axis, which 
thus is the strained c axis of the tetragonal lattice (Fig. [2]), may be answered as follows. DFT 
calculations show that there is a strong energetical preference for two nearby B atoms to oc¬ 
cupy the same type of octahedral interstitial position. The substrate may act as a seed layer 
to establish that the Erst B atoms occupy octahedral interstitial in direction perpendicular 
to the substrate. The DFT results then suggest that such domain will grow preferentially, 
compared to a random interstitial occupation or formation of multiple domains, both hav¬ 
ing higher energy. Multiple variants with differently aligned c axes are unfavorable due to 
the required interface, i. e. twin boundary, energy. Thus, Fe-Co-B Elms exhibit a quasi 
single crystalline state with one preferentially oriented strained c axis. This axis aligns 
perpendicularly to the Elm surface due to the square surface symmetry of the Au-Cu bufier. 

As second main result, the structural measurements revealed spontaneously strained 
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Feo. 38 Coo. 62 -B phases similar to those in Fe-Co-C^. The 20 nm thick films are much thicker 
than coherently strained films. The tetragonal strain is not dependent on film thickness: 
100 nm thick dims with the same compositions (see supplementary^) exhibit the same c/a 
ratios as the films presented here. The studied spontaneously strained Fe-Co-B films con- 
hrm our DFT calculations, which suggested minima of total energy depending on the B 
content iTable [IT|) . However, and this is an important difference to the Fe-Co-C films, the 
spontaneous strain indeed depends on the particular B content - at least for low B contents 
up to about 4 at%. 


B. Magnetic properties of the Feo.38Coo.62-B films 

In order to compare the magnetic properties of spontaneously strained Fe-Co-B films, the 
anisotropy constant iFu of the strain related uniaxial MCA was determined from hysteresis 
measurements as described in Ref.-. The results are summarized in Figure [6](b). Similar 
to c/a, K\] also has a maximum in dependence of the B content. The highest MCA with 
K\] = 0.53 MJ/m^ is reached in the 2.5 at% B containing Feo. 38 Coo .62 film. Films with smaller 
B content have a lower MCA because of their reduced tetragonal strain. B contents above 
2.5 at% in Feo. 38 Coo .62 do not increase the MCA, although the c/a ratio’s maximum is at a 
B content of 4.2 at% (Figure [6](a)). The observed decrease of MCA with further increased B 
content occurs simultaneously with a significant reduction of saturation magnetization in the 
films (Figure [T]). While our theoretical results described a reduction of Ms by about 1.5% 
with each percent B (open stars), the experiments (squares) reveal a stronger dependence 
between B content and magnetic saturation; The decrease is about 2.5 % with each 1 at% 
B. A possible reason might be a decreased Curie temperature due to the addition of B. 
However, this is not expected^ for B additions in Fe. We rather attribute the higher 
dependence of Ms on the B content to a stronger perturbation of the Fe-Co lattice by the 
B atoms in the PLD prepared films than it is described by the CPA approach, where the B 
atoms occupy regularly distributed interstitials. The calculations already revealed that Fe 
and Co atoms have a strongly reduced magnetic moment, when B is their nearest neighbor. 
In comparison to our experimental study, magnetron sputtered Fe-Co-B films^, where the 
material was deposited with substantially lower kinetic energy, exhibit a lower decrease of 
Ms with increasing B content, namely 1.5% per added 1 at% B, which is consistent with 
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FIG. 7. (Color online) Magnetic saturation in dependence of the B content for the 

Feo. 38 Coo. 62 -B films (squares) and the corresponding calculated structures (stars) according to 
Table |TI1 Linear fits have been added as broken lines. 


our theoretical results (Table HTj) . 

We argue that both, the strong reduction of Ms and the decrease of MCA observed in our 
samples originate from an interplay of the proceeding amorphization due to the increasing 
B content and the high energy impact of the ions during PLD^^ which may disturb the 
local environment of the Fe and Co atoms drastically. However, when we compare the 
theoretical and the experimental results of K\j (Figure [H](b)), we find increased values for 
the Feo. 38 Coo. 62 -B thin dims at low B contents. We attribute this discrepancy to a lack of full 
potential effects in DFT. At B contents above 4at%, this underestimation is overcome by the 
discussed structural effects. We thus observe a very good agreement of the experimentally 
measured magnetic anisotropy with the theoretically determined value for (Feo. 4 Coo .6 )24B, 
but no longer for (Feo. 4 Coo. 6 )i 6 B, i. e. 6at% B (open circles in Figure |6](b)). 

Concluding the magnetic measurements of the Feo. 38 Coo. 62 -B films, we observe a direct 
dependence of the magnetocrystalline anisotropy on the spontaneous strain for low B con¬ 
tents. However, an increased tetragonallity in the Fe-Co-B lattice did not necessarily imply 
an increased MCA for the complete films due to the formation of non-crystalline or disturbed 
regions. 
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FIG. 8. (Color online) Measured tetragonal lattice strain (full squares) and uniaxial magnetocrys¬ 
talline anisotropy constants K\j (full circles) in dependence on Fe-Co ratio in 20 nm films with 
4at% B on Auo, 55 Cuo, 45 . DFT results for c/a and K\j are added (open stars). 

C. Structural and magnetic properties in dependence on the Fe/Co ratio 

The study of various Fe/Co ratios by DFT calculations (Table HB motivated for a com¬ 
parison of films with unaltered B content. We chose Fe^jCoi-^; dims with 4at% B due to 
the observed strain maximum (Figure [6](a)). The measured tetragonal strain and the ob¬ 
served perpendicular anisotropy constants K\] in these 20 nm thick films are summarized in 
Figure [HI With regard to the tetragonal strain, we observe a slight decrease with increasing 
Fe content. Although there is only little variation between the films, this hnding is qualita¬ 
tively consistent with the theoretical predictions for (Fea,Coi_a,) 24 B in Table HIl which also 
give higher c/a ratios for higher Co contents. Due to already described reasons, we do not 
reach as high strains as calculated by DFT. The highest experimentally observed c/a is 1.045 
for X = 0.36. 

When comparing the magnetocrystalline anisotropy energies K\], we do not observe a 
strong dependence on the Fe/Co ratio, which in principle fits to the weak dependency 
proposed by DFT (see Figure [8]). K\j is about 0.5MJ/m^ on average, but slightly lower 
on the Co rich side and slightly higher in the dims with more Fe than Co. The maximum 
observed K\j is 0.54 MJ/m^ for x = 0.6. The DFT calculations, however, predicted the 
highest MCA for the lowest Fe content. 

From the Fe^jCoi-^, dims with 4 at% B we thus conclude that higher Fe contents {x = 0.6) 
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not only lead to a higher magnetic saturation Ms, but are also benehcial for a high MCA 
at low tetragonal strain. Such Fe rich films or additional concepts to stabilize higher strains 
beyond c/a = 1.05 obtained here are considered promising for further research. As one way 
to fulhll the latter task, we suggest the choice of a preparation method with less energetic 
impact, where the effect of lattice perturbation due to the B atoms might be smaller. 


V. CONCLUSION 

In this combined theory-experimental study, Fe-Co-B is introduced as an alloy with spon¬ 
taneous strain and hard magnetic properties. Structure calculations based on EMTO predict 
strong tetragonal lattice distortions up to c/a = 1.25 in dependence on the B content, if all 
B atoms occupy octahedral interstitials along the c axis. This arrangement is preferential, 
when compared to an occupation of other possible interstitial sites. SPR-KKR results for 
the distorted Fe-Co-B lattices give strong magnetocrystalline anisotropies when the former 
cubic symmetry is broken. For llat% B, K\] should be comparable to shape anisotropy. 
Experiments, however, show that the amount of B atoms which strains the Fe-Co lattice is 
limited. At a B content higher than 4at%, the strain decreases. A supersaturation of the 
octahedral lattice sites with B atoms is considered to be responsible for this deviation from 
the ideal behavior of B interstitials as treated in DFT. For low B contents, our theoretical 
results are well confirmed by experiments, in particular when regarding K\j. The highest 
magnetocrystalline anisotropies are observed for Feo. 38 Coo .62 with 2.5 at% B or in Fe richer 
Feo. 6 Coo .4 with 4at% B. K\] of these alloys is above 0.5MJ/m^. Although this value does 
not exceed the shape anisotropy, the spontaneous strain should allow for a preparation of 
Fe-Co-B with an easy axis of magnetization, which is dehned by structure and not by shape. 
The approach to exploit spontaneous strain and the related MCA thus gives much more 
opportunities than the substrate induced strain in Fe-Co, where the strain is limited to 
ultrathin films. 
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